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The properties of Josephson devices are strongly affected by geometrical effects such as those asso- 
ciated with the magnetic field induced by the bias current. The generally adopted analysis of Owen 
and Scalapino [Phys. Rev. 164, 538 (1967)] for the critical current, /c, of an in-line Josephson tunnel 
junction in presence of an in-plane external magnetic field, He, is revisited and extended to junctions 
whose electrodes can be thin and of different materials. We demonstrate that the asymmetry of 
the magnetic diffraction pattern, /c(i/e), is ascribed to the different electrode inductances for which 
we provide empirical expressions. We also generalize the modeling to the window-type junctions 
used nowadays and discuss how to take advantage of the asymmetric behavior in the realization of 
some superconducting devices. Further we report a systematic investigation of the diffraction pat- 
terns of in-line window- type junctions having a number of diverse geometrical configurations and 
made of dissimilar materials. The experimental results are found in agreement with the predictions 
and clearly demonstrate that the pattern asymmetry increases with the difference in the electrode 
inductances. 

PACS numbers: 74.50. +r,85.25.Cp,74.25.N- 
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I. INTRODUCTION 

In the first years after the discovery of the Josephson 
effect it was realized that, as the size or the critical cur- 
rent of a planar Josephson tunnel junction increases, the 
influence of the magnetic field induced by the Joseph- 
son current itself becomes more and more importaniPS^ 
Significant advances in the modeling of the so-called self- 
field effects were made by Owen and ScalapincP (OS) 
who considered the geometrical configuration most suit- 
able to analyze the self- field, namely, the one-dimensional 
in-line geometry shown in Fig. [TJa) where the externally 
applied bias current, /, flows in the direction of the long 
dimension, L. The figure also shows the coordinate sys- 
tem used in this work. Long Josephson tunnel junctions 
(LJTJs), i.e., junctions whose dimension L perpendicular 
to an externally applied magnetic field, i^e, is large com- 
pared to the Josephson penetration depth, Aj, behave 
like an extreme type-II superconductor; they exhibit a 
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Meissner effect in weak magnetic fields, and vortex pene- 
tration starts at a critical field He. This behavior can be 
probed into by studying the magnetic field dependence 
of the maximum tunneling supercurrent, /c, because the 
Meissner region and the vortex structure are reflected 
in the Ic versus He curve in a very characteristic way. 
The strength of the current-induced field is measured 
by the ratio of L/Aj. The OS analysis was restricted 
to the ideal case in which the junction electrodes had 
the same width and were made of the same bulky ma- 
terial. Despite these severe restrictions the OS modeling 
has been used without further consideration the design 
of any device based on LJTJs. One of the purpose of the 
present work is to generalized the OS theory taking into 
account geometrical or electrical differences in the elec- 
trodes. We will extend our interest also to junctions fab- 
ricated with the so-called whole wafer processes, which 
make use of a tri-layer structure from which individual 
junctions are later defined. Due to the large diffusion of 
these reliable processes, window junctions have become 
widespread while, in contrast, naked (or bare) junctions 
are a rarity. While the interaction of the Josephson tun- 
nel junction with its embedding circuitry has received an 
exhaustive and adequate attentionP®, the changes due 
to the self- field effects still remain an unexplored topic, 
albeit the exact knowledge of the current induced field is 
highly desirable for the realization of superconductor de- 
vices based on LJTJs such as, for example, oscillator^, 
magnetic sensors^^ and rectifiers-*^^. Besides modeling, 
a thorough experimental investigation of the self-field 
effects has been carried out for window-type A/'6-based 
LJTJs having in-line configuration and electrodes of dif- 
ferent widths, thicknesses and materials. The results are 
well aligned with the expectations and demonstrate that 
window-type LJTJs can be designed with a greater flex- 
ibility in the control of the self-field effects. The issue of 
the asymmetry in the threshold curves of Josephson de- 
vices was introduced long ago and its importance was rec- 
ognized in determining the performances of amplifier^^^ 
and magnetometers^^. Currently, the search for new and 
better ways t o achieve and control the asymmetry is still 
on-goin^Sf^. 



II. MODELING 

A. Background review 

In its simplest form, the area of a planar Josephson 
tunnel junction is defined by the overlap of two supercon- 
ducting films with rectangular cross section and weakly 
coupled through a thin tunneling barrier. Fig. [TJb) de- 
picts the vertical cross section of an in-line Josephson 
tunnel junction. An in-plane external magnetic field, H^^ 
is applied in the F-direction, i.e., perpendicular to the 
junction length. The very thin insulating layer between 
the superconducting films has length L and width W (not 
shown). Ih{X) and It{X) denote the local supercurrents 
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FIG. 1. (a) Schematic of an in-line junction with symmetric 
current bias, /, in the presence of an external magnetic field, 
ife, applied along the y-axis.; (b) Simplified vertical cross 
section (drawn not to scale) of a in-line Josephson tunnel 
junction. The junction top electrode is in gray, while the base 
electrode is in black; the tunneling insulating layer in between 
is hatched. Also shown is the coordinate system used in this 
work. 



flowing, respectively, in the bottom and top junction elec- 
trodes within a distance of the order of the penetration 
depth from the film surfaces and parallel to the insulating 
layer. X G [— L/2,L/2] is the laboratory spatial coordi- 
nate. Throughout the paper the subscripts h and t refer 
to the base and top electrode, respectively. Further the 
currents are positive when they flow from the left to the 
right. In the vast majority of practical cases the junction 
electrodes are comparable or thinner than their penetra- 
tion depths and the currents can be well assumed to be 
uniformly distributed over the film cross section. and 
It also take into account the screening currents, /5c, that 
circulate to expel the magnetic field from the interior of 
the junction. 

For in-line LJTJs, it is important to distinguish between 
the symmetric and asymmetric biasing: in the former, 
the bias curren t, /, enters at one extremity and exits at 
the otheiEiniBt /,(L/2) = h{-l./2) = I and /t(-L/2) = 
/5(L/2) = 0, while in the latter, t he bias c urrent enters 
and exits from the same extremit}fflIlIinilSl /^(— L/2) = 
-h{-L/2) = I and h{L/2) = /,(L/2) = 0. 

The gauge- invariant phase difference (j) of the order pa- 
rameters of the superconductors on each side of the tun- 
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nel barrier obeys the Josephson equation^^ 

Jz{X) =J^smcf>{X), 
kV(/)(X) =Hxri, 



(la) 
(lb) 



in which Jz is the Josephson current density and Jc is 
its maximum (or critical) value, which is assumed to be 
uniform over the barrier area. Eq.(lb) states that the 



phase gradient is everywhere proportional to the local 
magnetic field H and parallel to the barrier plane. Here 
n is the versor normal to the insulating barrier separat- 
ing the superconducting electrodes and n = <l>o/27r/io<^m 
has the dimension of a current (^o is the magnetic fiux 
quantum and /io is the vacuum permeability) . Neglecting 
the insulating barrier thickness^^. 



dn 



Xh tanh 



db 
2Xh 



Xt tanh 



dt 
2Xt 



(2) 



is the junction magnetic thickness in and A^^^ and d^^t are, 
respectively, the bulk magnetic penetration depths and 
thicknesses of the films. Eq.([2| reduces to dm = X^ -\- Xt 
in the case of thick superconducting films (^5^^ > 4A5^t). 
The net current, /, crossing the tunnel barrier is given 

by 



/•L/2 

I = W Jz{X)dX. 

J -L/2 



(3) 



Throughout this paper we will limit our interest to the 
zero- voltage time-independent state; this can be achieved 
as far as / is smaller than the junction critical cur- 
rent, /c. We also assume that the junction wid th, W, 
is small er than the J osephson penetration deptli^^^^^^, 
A J = ^y^o/27^^^odjJc^ setting the length unit of the phys- 
ical processes occurring in the Josephson junction; here 
dj is the junction current thickness^^ (see later) 



d. 



o — A5 coth ^ 
A 



db 

b 



Xt coth ^ > dr, 
Xt 



(4) 



For thick film junctions, dj = dm = Xi^ -\- Xt. By its 
definition, = {dj/dm)Jc^'j- 



It is well known ^ ^ that combining Eqs.(la) and (b) 
with the static Maxwell's equations, a static sine-Gordon 
equation is obtained that describe the behavior of a one- 
dimensional in-line LJTJ 



A 



2^ 



:sin0(X). 



(5) 



Equation([5| was first introduced by Ferrel and Prange^ 
in 1963 in the analysis of asymmetrically biased in-line 
LJTJs; few years later, reported an extensive study 
of its analytical solutions in terms of elliptic functions 
for symmetrically biased in-line junctions (provided that 



L > ttXj/2). As reported by several author^I^SlMl^ the 
largest supercurrent carried by a very long in-line LJTJ 
is 4/0, where /q = JcWXj is a characteristic junction cur- 
rent (generally from a fraction of a milliampere to a few 
milliamperes) and depends on the junction normalized 
length, £ = L/Aj, as /o(f) = /otanh^/2. 
As already stated, for small fields a LJTJ behaves as a 
perfect diamagnet by establishing circulating screening 
currents which maintain the interior field at zero. This 
Meissner regime is refiected by a linear decrease of Ic 
with He. The threshold curves, Ic{He)^ have been the 
subject of many analytical, numerical and experimen- 
tal work^^ M^. It is well known that the critical mag- 
netic field i^^Hc = ^o/ir/ioXjdm = 2JcXjdj/dm- The 
gradual crossover of He from short (-^ << 27r) to long 
(^ > 27r) junctions has been numerically computed^^' 
and it was found to be well described by the empirical 
relationship: Hc{£) = He cothi/ir in very good agree- 
ment with the experimental findings reported in Ref.l^. 
Note that in the small junction limit, i 0, we re- 
cover He = ^0/ jiioLdm and the Fraunhofer-like magnetic 
diffraction pattern (MDP). 



B. Boundary conditions 

In this section we will derive the boundary conditions 
for an in-line junction. Let us define the sheet induc- 
tances: 



b,t 



_ MoA6,t 


tanh 


db,t 


2 


2Xb^t 


_ MoA6,t 


coth 


db,t 


2 


2Xij^t 



(6a) 
(6b) 



and observe that for any value of the ratio da/Xa^ it is 
Ca > fioXa/2 > Ta with a = b^t; in the thin film limit 
Ta = jJ^oda/^ and Ca = ji^X^/da is the film kineticP^ 
sheet inductance due to the inertial mass of mobile charge 
carriers (Cooper pairs). Furthermore^^, Ca ^ Ta = 
/ioAaCoth da/Xa, while Ca-Ta = liioXaCschda/Xa. Then, 
according to Eqs.([2| and Q, /nodm = 2(T5 + Tt) and 
jj^odj =Tb^Tt^Cb^ Ct. 

In the case of a Josephson structure made by two films 
of the same width, Weihnact^^ provided a general 
expression for the phase derivative as the sum of three, 
in general, independent terms 



<I>o d4> 
2^dX 



= jl^dmHe 



W 



{Ct^n)-^{C,^Tt). (7) 



More precisely, if we introduce the magnetic fiux changes, 
A<I>j, in the Josephson barrier proportional to the 
changes of the Josephson phas^f, A<l>j = <l>oA^/27r, 
each term in Eq.Q represents a magnetic fiux per unit 
length along the X-direction. Let us assume first that 
the junction is unbiased, so that the screening current. 
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/gc, only circulate in the electrodes to expel any exter- 
nal field, i^e, from the interior of the junction; charge 
conservation requires It{X) = —I}){X) = Isc{X). Then: 



^0 



(8) 



By resorting to the definition of inductance as fiux per 
unit current!^, we can identify Cjtl = f^odj/W as the 
inductance per unit length for the screening current^^, 
where the subsript JTL stands for Josephson Transmis- 
sion Line. The same expressioiP^"^ was obtained for a 
strip- type superconducting transmission line (STL) con- 
sisting of a dielectric layer sandwiched between two su- 
perconductors of width W; here, the current fed into the 
top electrode returns back in the base electrode acting as 
a ground plane, namely, hi^) = —It{X). In such a way, 
all fields are identically zero outside the waveguide (the 
only difference being that, for JTLs, the dielectric thick- 
ness can be neglected). The screening current in Eq.([8| is 
negative for positive external field, and vice versa. Since 
jioHe is the externally applied magnetic fiux density, B^^ 
and iiqIsc/W is the magnetic fiux density, Bsc^ induced 
in the barrier by the current I so Eq.([8| help us to un- 
derstand why we named dj^ and dj as the junction mag- 
netic and current thickness, respectively. Likewise /io<^m 
and fi^dj can be seen as, respectively, the magnetic and 
current sheet inductances of the junction. To find the 
boundary conditions for Eq.([5]) when the in-line junction 
is biased is not an easy task, because it requires the sep- 
arate knowledge of the properties of each electrode. All 
previous analytical approaches'^ dealt with transmis- 
sion line structures in which the current fiowing into the 
top electrode returns back in the lower electrode acting 
as a ground plane; in such cases the inductance per unit 
length of the transmission line is the only required pa- 
rameter. However, for our purposes we have to find the 
inductances, £5 and Ct per unit length of of the bottom 
and top electrodes, respectively. Resorting to Eq.Q, we 
easily identify them as: 



w 



(9) 



so that C})^ Ct = Cjtl'- this is a well known identity in 
the theory of two-conductor transmission lines^^'^^'^. If 
both films are thick, Ci) = Ct = /J^o{Xb + ^t)/'^W; in the 
opposite case, Cb^t = /io^^ ^/Wd^^^, that is when both 
films are thin their inductance is essentially of kinetic 
origin^. Classically, the bottom and top inductances 
have been merged in their parallel combinatioiP^'^ so 
that the role played by each supercurrent separately 
was lost; however, to correctly describe a biased in-line 
Josephson junction it is mandatory to keep the distinc- 
tion, since, each electrode transports its own supercur- 
rent, and, in general, Ib{X) ^ —It{X). The first term 
in Eqs.(|9| takes into account the magnetic and kinetic 
energy in the film carrying the current, while the second 



term is the magnetic energy stored in the opposite elec- 
trode. For a thin-film (symmetric) transmission line the 
ratio of the internal magnetic to the kinetic energy goeJ^ 
as (2(i/A)^/12, i.e., for thin films the kinetic energy dom- 
inates. At this stage we are still allowed to neglect the 
magnetic energy stored in the very thin insulating tun- 
nel barrier. It is worth to stress that Eqs.(|9| are valid 
as far as the energy of the magnetic fields outside the 
transmission line is negligible. 

The case of a superconducting strip of width W carrying 
a current / over a a superconducting shield was consid- 
ered a long ago'^. By resorting to the theorem of images, 
it was shown that, if the strip is close to a thick supercon- 
ducting shield, the field between them is approximately 
uniform and equals to H = I /W and, outside the edge 
of the film, it falls to zero within a distance of the or- 
der of the strip-to-shield distance. The shield screen- 
ing current, whose integral is equal to /, must therefore 
be uniformly distributed over that portion of the shield 
surface covered by the film. On the other side of the 
current-carrying strip the field is zero. In the absence of 
a shield, it is easy to show that the field on both sides 
of a freestanding strip will be equal and opposite and of 
mean magnitude H = I/2W. Hence, by bringing a su- 
perconducting plane close to a current-carrying strip, we 
double the field between the plane and the strip and re- 
duce it everywhere else. This effect can be used to reduce 
the effective inductance of superconducting elements by 
depositing them on top of another insulated supercon- 
ductor. These considerations allow the use of Eqs.([9| 
when a current locally fiows either in the bottom or in 
the top junction electrode, as it occurs in a symmetri- 
cally biased in-line junction. However, if in a given place 
Ib{X) and It{X) are comparable and fiow in the same 
direction, then a considerable magnetic field is induced 
in the outer space (unless the films are thin). 



C. Symmetric biasing 

To obtain the local magnetic field, i7y(X), we have to 
divide Eq.Q by ji^dm] in particular, at the boundaries 
of a symmetrically biased LJTJ we have 



dcj) 



dX 



-.HY[--\=He--^ = He-^-^ — , 
2 / l^odm dm ^JTL 



Ctl 

l^odm 



--H, 



dj 



For Ch = Ct 



Cjtl/'^ (and dj 



dr, 



im Cjtl W ' 
(10) 

), we recover the 



symmetric OS boundary conditions^ that were generally 
adopted ther eafte r for untrue symmetry reasons. In the 
early eightieP^^, the reported asymmetric behavior of 
samples that were believed to be symmetric led many 
experimentalists to abandon the in-line geometry in favor 
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of the overlap one. Indeed, the last term in Eqs. (10) is 



the magnetic field induced in the tunnel barrier by the 
current / fiowing in the lower or in the upper electrodes. 
Ampere's law applied along the barrier perimeter in the 
X-Y plane requires that the magnetic fields at the two 
ends of the junctions differ by the amount of the enclosed 
current: / = W^Hyil^/T)- i^W-L/2)]. It is easy to 
show that Eqs.([la|), ([3|, ^ and (p^ fulfih Ampere's law. 
From the boundary conditions (|10[), it follows that 



2H, 



dj I Ct 



dm W CjTL 



The last term, vanishing when Ch = Ct^ has been omit- 
ted in all previous analysis of LJTJs. The difference, 
Ct — Chi ill the self-inductances of the upper and lower 
junction electrodes is responsible of the distortion of 
the experimental curve^; it equals /io[AtCsch (d^/At) — 
Afocsch (d^/A^b)] and reduces to ii(^{\t — A5) in case of 
thick electrodes. We identify the dimensionless parame- 
ter a = {Ct — Cb)/CjTL as a direct measure of the asym- 
metry of the system, with — 1 < a < 1. The asymmetry 
can be ascribed to differences in the electrode thicknesses 
and/or materials, that is, it can have geometrical and/or 
electric origins. 



We like to point out that Eqs. (10) are very general and 

'ef- 



can be used to correctly describe the so-called se{ 
fects occurring in any substantially in-line LJTJ in which 
the bias current or a part of the bias current enters or 
leaves the junction from one or even both its ends. Un- 
fortunately, their implementation requires the separate 
knowledge of the bottom and top electrode inductances 
per unit length (rather than just their sum). 



the above reasons, in the rest of this paper we will only 
consider the more interesting case of the symmetrically 
biased in-line junction with no ground plane. 



E. Normalized units 

In normalized units of x = X/Aj, the differential equa- 
tion ([5| becomes: 



(l)a:x = Sin^(x), 



(12) 



with X G [— ^/2,^/2]. Normalizing the currents to Iq = 
JcWXj and the magnetic fields to Hc/2 = JcXjdj /dm^ 
the boundary conditions (10) for a symmetrically biased 
LJTJ read 



hi he. 



Chi 



Cti 



(13) 

With these notations, the normalized critical magnetic 
field he of a short Josephson junction is 2i\ ji. 



1 . Magnetic diffraction pattern 

Setting hi and hr at their extreme values ±2 in 
Eqs. (13), we obtain the normalized MDP, ic{he)^ in the 
Meissner regime 



D. Asymmetric biasing 

In the case of asymmetric bias, the boundary condi- 
tions are 



2 + 



Hy 



dm CjTL W 



(11) 



and, recalling that Cb -\- Ct = Cjtl^ we end up with the 
boundary conditions first found by Ferrel and Prange^ for 
thick electrode LJTJs. We observe that the difference in 
the inductances does no t play any role in this case. Many 
pioneeristic experiment j-*^^ * ^^^ were also carried out with 
in-line junctions built over a large and thick supercon- 
ducting ground plane that makes the current to fiow 
in the lower skin layer of the bottom electrode provided 
it is a thick film and is shielded by the upper part. When 
this is the case, £5 must be set to zero in Eqs. (10) or (11), 



z{he) 



for -2 < he < hmax 
for hmax <K <2, 



(14) 



with hmax = 2a being the field value which yields the 
maximum critical current ic{hmax) = 4. The second-last 
equality turns out to be very useful in the experiments to 
determine the asymmetry parameter from the analysis of 
the junction MDP: 



(15) 



SO that the biasing configuration becomes irrelevant. For 



Fig. [2] shows the theoretical Ic{He) for a very long in- 
line Josephson junction in the Meissner regime and for 
different values of a; the junction critical current, /c, is 
normalized to /q and the critical magnetic field, is 
the theoretical field value that fully suppresses the critical 
current. We remark that the patterns are antisymmetric, 
Ic{—He) = —Ic{He) and piece wise linear. They have 
different absolute slopes \dIc/dHe\ on the left and right 
branches, respectively, equal to jHodm/Cb and /io<^m/A- 



6 



4 
3 
2 

c 


-1 
-2 
-3 
-4 













/ 


















\ 










































\ 










































s 














/ 


























\ 












/ 

/ 














a 














\ 






/ 


/ 


























> 










/ 


/ y 
















-— -0.8 
- - -0.4 



-- 0.4 








































































































































08 






























































\ 


























/ 
















\ 






























































V- 




















\ 


















/ 













-1.0 



-0.5 



»"h/h"' 

e c 



1.0 



FIG. 2. Theoretical magnetic diffraction patterns Ic{He) 
of a very long in-line symmetrically current-biased Joseph- 
son junction in the Meissner regime for different values of 
At = Ct/CjTL- The critical current, /c, is normalized to 
^0 = JcWXj and the externally applied magnetic field, He, 
to the critical magnetic field He. 



2. Current diffraction 'pattern 

Let us consider now the case when the junction criti- 
cal current is modulated by the transverse field induced 
by a stationary current entirely flowing in either the top 
or bottom electrode. This so-called control current tech- 
nique has been used to produce local magnetic fields for 
digital applications of Josephson circuits since 196&^^. If 
the control current, /, is injected into, say, the top elec- 
trode, then the magnetic field at the inner surfaces is 
Ht = Cti / iJi^dmi as in Eq.([lQ|. The value of the control 
current for which the junction critical current Ic vanishes 
is Ic = ^o/nXjCt and will be named the critical control 
current. Likewise, Imax will be that value of / which 
maximizes the critical current, /c- In normalized units, 
ht = 2Ht/Hc = Cti/ CjTLi where i = I /Iq. Replacing 



with ht^ Eqs.(14) become 



^(^c + for —ic < i < in 



for L 



<i<ic, 



where tc h/h 2CjtlI C.t and l^ax ^max 
f^d^b — ^t)/^JTL- It is worth mentioning that 



(16) 



//o = 



-^max -^maa; 



(17) 



indicating that the degree of asymmetry is the same for 
the magnetic and the current diffraction pattern (GDP). 
In other words, as far as the modulation of the critical 
current concerns, a control current is equivalent to an 
externally applied field. It is then possible to identify the 



GDP slope with a current gaiiP^^ g = dlc/dl = dic/dt. 
From Eq.(IT7l), it follows 




X, ^ i ^ ic 



(18) 



With Ct > then t^ax is negative and so, in zero 
external field, we have a unitary current gain. However, 
large current gains can be achieved with samples having 
a large Ct/^b ratio by control current biasing the upper 
electrode to have t < tmnT- 



III. WINDOW-TYPE JUNCTIONS 

Weihnact's Eq.Q was derived considering that the 
thickness of the tunnel barrier is much smaller that typ- 
ical penetration depths so that the magnetic field in be- 
tween the electrodes is in the y-direction, is uniform 
along Z, and only depends on the X coordinate. Ac- 
cordingly, there are no fringing effects and, at the same 
time, the magnetic energy stored in the dielectric layer 
can be neglected. In addition, it was assumed that the 
base and top electrode match the width of the barrier, 
Wb = Wt = W. However, these conditions are not 
fulfilled in nowadays window- type planar tunnel junc- 
tions whose electrodes have quite different widths; typi- 
cally, Wi) > Wt > that is, two strips, of total width 
Wi = Wt — W^ exist along the junction sides where the 
top electrode overhangs the base electrode, but no tun- 
neling is possible due to the thick insulating layer. In 
case Wt > Wb, then Wi = W^ — W; however, through- 
out this paper we will assume that the base electrode 
is wider or at most matches the top electrode. The 
so-called idle region is formed after the patterning of 
the wiring film which provides the electrical connection 
to the junction top electrode. In this region surround- 
ing the tunnel area the insulation between the bottom 
and top electrode is provided by an oxide layer typi- 
cally made of a native anodic oxide of the base electrode 
and/or a deposited SiOx layer. The total thickness of 
this passive layer, dox^ is comparable or even larger than 
the electrode penetration depths, A^^t, and might also 
be comparable with the junction width, W. Then the 
screening currents in the upper and lower electrodes dis- 
tribute over an effective width larger than the junction 
width. This leads to a new smaller junction inductance 
per unit length C jtl corresponding to the so-callecPM^ 
inflation of the Josephson penetration depth occurring 
in window-type Josephson tunnel junctions. In Refs."^, 
each electrode was modeled as a parallel combination of 
two stripes having quite different oxide thicknesses re- 
sulting in a rather involved expression of the effective 
current thickness. From a phenomenological point of 
view, a unidimensional junction with a lateral idle re- 
gion behaves as a bare junction having the same width, 
W ^ and length, L, but with an effective current thick- 
ness, d'j^ given by a proper combination of the current 
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thicknesses of the naked junction, dj and of the idle re- 
gion, di, namely: dj = dj/[l + {dj /di){Wi/Wj)] < dj, so 

that C jTL 



IJLod'j/W and \'j = ^^of^TTfiod'jJc > Xj. 

The prime symbol (') labels the parameters relative to 
the window junction. Taking into account the thickness, 
d^, of the wiring layer, it is: 



dj — Xb coth ^ 
Xb 



At coth 



{dt ^dyj) 

Xt ' 



dj = Xh coth ■ 



A. 



Xt coth 



A, 



Outside the junction area the idle region takes the 
form of a microstrip-line made by two electrodes of fi- 
nite width and thickness. Chang^^ considered the case 
of a superconducting strip transmission line, i.e., a struc- 
ture consisting of a finite-width superconducting film of 
thickness h over an infinite (and thick) superconducting 
ground plane. As far as the strip linewidth, exceeds 
about the insulation thickness, t, the inductance per unit 
length, jCstl^ was analytically derived as 



-STL 



W t 
h'h 



with the fringing- field functional first introduced in 
Ref.!^, being always larger than unity; the fringing fields 
have the effect to increase the system energy, i.e., to re- 
duce the inductances per unit length. K decreases with 
the the ratio w/h and increases with t/h. Definitely, 
Chang's results can be used when Wb » Wt, but, unfor- 
tunately, no analytical expression is available for a STL 
when both electrodes have finite and comparable widths. 
The situation complicates if Ib{X) 7^ —It{X). Therefore, 
we resort again to Eq.([7| to determine the flux per unit 
length just outside the tunneling area 

= fiod'r^H, + C'tit - C',h. (19) 

In other words, for a window- type LJTJ the left (right) 
boundary condition is determined by the property of 
the longitudinal idle region to the left (right) of the left 
(right) junction end. The longitudinal idle region also 
acts on the junction as a capacitive load, which does not 
play any role, as far as the static properties concern. We 
propose that, provided that Wb slightly exceeds Wt^ the 
inductances per unit length of the bottom and top elec- 
trodes at the extremities of a symmetrically biased in- 
line junction (where the bias current, /, at the junction 
extremities either flows is the bottom or in the top elec- 
trode) are 



ox 

Cb-\-Tt -\- fiodox 
KbWb ' 



(20a) 
(20b) 



where Kb^t are fringing field factors that remain to be 
determined: with Wb > Wt we expect Kb > Kt ^ I. 
Eqs.( |20a] ) and (b) rely on the fact that a current It {h) 
flowing in the top (bottom) electrode induces a magnetic 
field It/Wt (h/Wb) between the electrodes. The last 
terms in each of the previous equations take into account 
the magnetic energy stored in the thick oxide layer. In 
Eq.(19) we also have considered that the magnetic thick- 
ness outside the tunneling area is = dm^dox (but this 
is a secondary field focusing effect^^^^). Consequently, 
for any window- type Josephson tunnel junction the crit- 
ical field is lower than that of a naked one; for a window- 
type LJTJ it is 



H'^ = ^o/iTfiod'^X'j < He. 



(21) 



Ultimately, the boundary conditions Eqs.(lO) still hold. 



if we replace d^, dj, Cb^t and Cjtl for the naked junc- 
tion with their respective counterparts, d^, d'j, C'b^t and 

JTL^ for the window junction. It should be clear that, 
this time, the inductance per unit length, C stl^ of the 
superconducting (non- Josephson) transmission line out- 
side the Josephson area is not given by the sum C'b^C't, 
otherwise the magnetic energy stored in the dielectric 
layer would be counted twice. Naturally, in general, 
C'b + C't not even matches C jtl- For the samples fab- 
ricated by means of a tri -laye r process, the top electrode 
to be considered in Eqs.( |20a'| ) and (b) is the wiring layer 
which does not necessarily cover the whole barrier area, 
so that it could as well be Wt < 

Indeed, the in-line approximation fails, if Wb or Wt 
are much larger than W, because part of the bias cur- 
rent enters the tunnel barrier also along the long junc- 
tion dimension, L, and a mixed in-line-overlap model^^'^ 
should be adopted in which the Josephson phase (j) obeys 
the time-independent perturbed sine-Gordon equation 



sin^(X) = A^ 



dX^ 



JcW 



Xov{X) is the distribution of the externally applied bias 

L /2 

current, Iqv = J_L/2-^ov{X)dX , giving the fraction of 
the bias current / that enters the junction through the 
long dimension (overlap component). Nevertheless, the 
junction configuration is still substantially in-line and, as 
such, the boundary conditions discussed so far still apply 
with the only caveat that the in-line component of the 
bias current is / — Iqv • The consequences resulting from 
the asymmetric boundary conditions imposed by a non- 
uniform external magnetic field at the extremities of both 
short and long Josephson junctions have been recently 
investigated^^ ^^; the field asymmetry is responsible for a 
degeneracy of the critical field, that was numerically 
demonstrated and experimentally verified. To better un- 
derstand this point, we sketch in Fig. [sj^a) the geometry 
of a 100 /im long and 1.2 /im wide, symmetrically biased 
window-type Nb-AlN-NbN junction with Wt = 4 jam 



(a) 




-500 500 

H(A/m) 
(b) 



FIG. 3. (a) Electrode configuration of a 100 /xm long and 
1.2 /xm wide, window-type Josephson tunnel junction hav- 
ing quite different electrode widths: Wt — 4 /xm and Wh — 
30 /xm; (b) Its experimental magnetic diffraction patterns, 
Ic{He), when symmetrically biased: Ic,max = 1.10mA, Hd = 
±580 A/m and Hc2 = ±830 A/m. 



and Wb = 30 /im. Fig. [sj^b) shows the MDP in which the 
two critical fields, Hd and Hc2^ are quite evident. Since 
in real measurement Ic never vanishes before the next 
lobe "grows up", the critical field values are obtained by 
extrapolating to zero the linear branches of the principal 
lobes; this is indicated by the gray dashed lines in this 
and forthcoming plots. A similar pattern asymmetry was 
first reported in Ref.!^ and was erroneously ascribed to a 
low uniformity of the barrier properties. However, this is 
an extreme case in which the asymmetry is mainly due to 
the very different widths of the electrodes. The experi- 
mental findings that will be presented in the next Section 
refer to window junctions having > 1.5Wt for which 
the degeneracy of the critical fields, although observable, 
is small enough to be ignored. 



/i; = K 



h' = K 



where we have introduced the reduced inductances A 



(22) 



6,t 



h,tl JTL and, in general, + AJ 7^ 1. Setting h\ and 
h!^ at their extreme values ±2 in Eqs.(22), we obtain the 
MDP i'dh^) in the Meissner regime lobe 



i'eiK) 



for -2<h',< h'^,, 
for h'^,, < K < 2, 



(23) 



with h'^^^ = 2{C\ - C't)/{C\ + C't) being the normal- 
ized field value yielding the maximum critical current 
i'cih'^a.) = 4/(A; + AJ) = ^C'jtlI{C\ + C\). Still we 
can define the asymmetry parameter 



a 



H: 



/ 

max 



JTL 



H' 



(24) 



and express the ratio of the inductances per unit length 
as 



C't 



1 



(25) 



when a' is positive then C t < CJ h and vice versa. If the 
modulating magnetic field is the transverse field induced 
by a control current, /, flowing in the top electrode, then 
R'l = C'tl/l^odm — d'jK[I / d'^W . In normalized units, 
h't = 2H't/H'^ = Kl', where t' = I/Tq is the reduced 
control current. Replacing h'^ with K[i' in Eq.(23), we 
get the CDP of a window- type junction 



where 



2/A; and 



for - 
for 



<l'<l' 



(26) 



1)/(AJ, + AJ). The current gain, g' 
window- type junction is 



h'^^JK = 2{A'JK - 
di'^/di' ^ for a in-line 



for —i 
for I' 



< < 



(27) 



Normalized units 



The boundary conditions for a window junction are ob- 
tained dividing Eq.(19) by jJ^od'^. In normalized units of 
x' = XI W J, with [-f /2,f /2] where f = L/A^^ is 

the normalized length of the window junction; normaliz- 
ing the currents to /q = Jc W^.i and the magnetic fields 
to E'J2 = JcX'jd'Jd'^, Eqs.([l3f become 



Since most real samples have Wh > Wt^ in the case of 
electrodes having the same penetration depth, then jC^t > 
£'5, suggesting that, for applications it is preferable to 
inject the control current into the top electrode which, 
typically, with respect to the bottom electrode, has a 
smaller width and so a larger inductance per unit length. 
Furthermore, can be increased by simply reducing the 
wiring width at the junction extremities. This can be an 
advantage for some applications. 
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IV. EXPERIMENTS 

The findings of the previous Section can be experi- 
mentally verified by measuring the magnetic and cur- 
rent diffraction patterns of symmetrically biased in-line 
window- type LJTJs. In particular, Eqs.(24) and (25) al- 



low the determination of the ratio of the inductances per 
unit length which can be compared with the expected 
value from Eqs.(20a) and (b): 



= a 



Wt Cb + Tt + Mocio 



(28) 



where a = i^^ /K^ is a factor which takes into account the 
film asymmetrjEH and is equal to 1 when the electrodes 
have the same geometrical and electrical parameters. We 
will test the validity of the Eq.(28), for samples having 
V^5 = Wt and = l.bWt. 



A. Samples 

In order to have a solid benchmark of data, we have 
compared the static properties of thin-film Nb-Nb and 
Nb-NbN samples having the same geometrical config- 
uration. The NbN films used in this study were de- 
posited by dc magnetron sputtering of a pure Nb tar- 
get in an argon and nitrogen sputtering gas. Although 
the samples were not heated intentionally, the substrate 
temperature slightly rose during deposition; neverthe- 
less surface temperature never exceeds 120 allowing 
for a lift-off process. The superconducting transition 
temperature (determined as the resistivity midpoint by 
a four-point method) for samples deposited under op- 
timal conditions is 15.5 K. Polycrystalline Niobium Ni- 
tride has a dirty-limit penetration depth"^^ j^^^^j^^" = 
4:.2K) = 370 nm, several times larger than that of epi- 
taxially grown Niobium, XNb{T = 4.2i^) = 90 nm. 
All samples were long window- type in-line junctions hav- 
ing physical length L = 100 /im and width W = 1.2 /im 
with a 1.4 /im wide idle region on each side, so that 
Wi = 2.8 fim. They were fabricated with the same pa- 
rameters for the deposition and anodization of the base 
electrodes, as well as for the deposition of the passive 
layer: db = 190 ± 10 nm and dox = 220 ± 10 nm. Further, 
the top and wiring layers had quite similar thicknesses: 
dt,Nb = dt^NbN = 65 ± 5nm, d^^Nb = 470 ± 10 nm and 
dw,NbN = 390 ± 10 nm. Therefore, the material used 
for the top (and wiring) electrode was the only notable 
difference. More importantly the wiring layers were Nb- 
Nb samples, duj^Nb ^ ^^Nb^ and thin for the Nb-NbN 
ones, d^^NbN ^ ^NbN^ in which case the inductance is 
predominantly kinetic. We also note that the thickness, 
dox^ of the passive layer is comparable to the electrode 
thicknesses. Table I reports the relevant electric param- 
eters (at 4.2 K) for the Nb/Alox/Nb and Nb/AlN/NbN 
Josephson tunnel junctions used for this work. For all 
samples it was L >> A- < W. All measurements were 



Sample Jc di dm dj Xj He d'm dj X'j 

kA/cm^ nm nm nm fim A/m nm nm /am A/m 
Nb-Nb 11 400 160 180 3.7 890 380 89 5.3 260 
Nb-NbN 3.6 780 270 530 3.6 520 470 210 5.9 190 



TABLE I. Relevant electric parameters (at 4.2 K) for the 
Nb/Alox/Nb and Nb/AlN/NbN Josephson tunnel junctions. 
The values are approximated to two significant digits and the 
prime symbol (') denotes the parameters relative to the win- 
dow junctions. All samples have the same length, L = 100 /xm, 
width, W = 1.2 /im, and idle region width, Wi — 2.8 /j.m. 



carried out at 4.2 K. The experimental set-up has already 
been described elsewhere^. 



B. Magnetic diffraction patterns 

The upper panel of Fig. |4] shows the geometry of a in- 
line junction realized by equal width electrodes, Wb = 
Wt = 7/im. The electrodes are vertically shifted, but 
they preserve the symmetry with respect to the junction 
axis. In Fig. we report the MDP of such junction 
when fabricated with the all-Nb technology. The small 
positive asymmetry, a' = H'^^^/H'^ ^ 7%, can be fully 
ascribed to the difference in the electrode thicknesses, 
d,„ ~ 2.5 db, and implies that C^t is sligthly smaller than 

0.90, 



C'b. From Eqs.(28) with a = 1, we have CJ t/ CJb 



in quite good agreement with the value of about 0.88 
obtained from Eq.(25). Fig. |4]^c) is the counterpart of 
Fig. Qb) for a Nb-NbN junction having the same ge- 
ometry. Now the asymmetry parameter, ol ^ —0.53, 
is negative indicating that C't > CJ b- The factor a in 
Eqs.pS]) must be set to 1.9 to reproduce the measured 
x^WoCjtlCJb = 3.3. 

Figs. [5|b) and (c) compare the MDPs of two symmet- 
rically biased junctions having the geometrical configu- 
ration, depicted in Fig. [sja) {Wb = 1.5Wt = 6/im), but 
made by different materials, respectively, Nb-Nb and Nb- 
NbN. We now have asymmetries of a' ~ —30% and 
—89%, respectively. For the former sample the asym- 
metry is mainly ascribed to the difference in the elec- 
trode widhts, while for the latter it is further enhanced 
by the large penetration depth of the upper (NbN) elec- 
trode. To reproduce these values, we had to set a in 
Eqs.([28| to 1.4 and 3.3, respectively. From Eqs.([25| we 
have C t,NbN ~ 4:.6C't,Nb ~ S.5C'b- 



C. Current diffraction patterns 

Fig. [6]^a) shows the CDP of the same Nb-NbN sam- 
ple of Fig. p[c) when the control current fiows in the top 
(wiring) electrode and in absence of an externally applied 
field. Comparing it to the MDP in Fig. [sjc), we observe 
a similar qualitative behavior and, as expected according 
to Eq.(17), the measured asymmetry parameters are the 
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FIG. 4. Experimental magnetic diffraction patterns, /c(ife), 
for two symmetrically biased in-line long Josephson tunnel 
junctions having equal width electrodes {Wb = Wt = 7/im), 
but made of different materials: (a) electrode configuration; 
(b) Nb-Nb sample: Ic,max — 2.28 mA, Hmax — ^OA/m and 
He = ±710 A/m and (c) Nb-NbN sample: Ic,max = 1.08 mA, 

ax 



-465 A/m and He = ±870 A/r 



< 

B 0-6 



0.4 
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FIG. 5. Experimental magnetic diffraction patterns, Ic{He), 
for two symmetrically biased in-line long Josephson tunnel 
junctions having the same geometrical configuration {Wt — 
Aiim and Wh = G/j^m), but made of different materials: 
(a) electrode configuration; (b) Nb-Nb; Ic,max = 2.23 mA, 
Hmax = -175 A/m and He = ±585 A/m and (c) Nb-NbN; 
Ic,max — 1.03 mA, Hmax — —5^5A/m and He — ±680 A/m. 



same within the experimental uncertainty of a few per- 
cent. Furthermore the current gains, measured by the 
slopes of the left and right branches of the Meissner lobe 
are in agreement with the expectation of Eq.(27). 



agree with Eq.(29). Other results in agreement with the 



If the control current is fed to the bottom, rather than the 
top electrode, the normalized inductance per unit length 
of the junction bottom electrode, AJ^, must replace in 
the numerator of the fractions in Eq.(26), so that 



for - 
for 



< ^" < 4 



(29) 



where now i'^ = 2/ A', and = 2(AJ/A; - 1)/(A; + AJ). 
Fig. erb) is the CDP of the same sample when the control 
current is injected into the base electrode. The asymme- 
try parameter is still unchanged and the pattern slopes 



theory (and not reported here) where obtained for the 
Nb-Nb sample of Fig. [sja). The wide ranece of linear- 
ity of the CDPs is very attractive for the realization of 
cryogenic current amplifiers with a large dynamic range 
especially because large slopes can be achieved in the 
[—f^c^f^max] interval. The asymmetry in the electrode in- 
ductance can help to significantly improve the gain of 
a current amplifier over a device with symmetric induc- 
tances. The most practical way to achieve this is to in- 
ject the signal current in the electrode having the largest 
inductance per unit length. In practice, the top one is 
chosen which can be made thin, narrow and of a ma- 
terial with a large penetration depth, such as polycrys- 
talline NbN; possibly, on the contrary, the base electrode 
should be thick, moderately wider and made of a material 
having a low penetration depth, such as epitaxial Nb. 
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FIG. 7. Schematic view of a long Josephson tunnel junction 
having a ring-shaped top electrode (in gray); the base elec- 
trode is in black and the tunnel area is white. 



Of course, the current-induced magnetic fields do not suf- 
fer from such screening. Therefore, a systematic inves- 
tigation of magnetic and current diffraction patterns of 
symmetrically biased in-line L JTJs could be useful to un- 
derstand the demagnetization effects in superconducting 
thin-film structures in presence of an in-plane external 
magnetic field. We observe that the inductance ratio 
(and so the asymmetry parameter) can also be extracted 
from the ratio of the top and bottom critical control cur- 
rents. 



FIG. 6. Experimental current diffraction patterns for the 
same sample of Fig.jsjjc). The control current / is injected: (a) 
into the top electrode: Ic,max = 1.03 mA, Imax = —0.46 mA 
and Ic = ±0.56 mA, (b) into the base electrode: Ic,max = 
1.03 mA, Imax = —3.4 mA and Ic = ±4.6 mA. g' is the cur- 
rent gain given in Eqs.(27) or (29). 



D. Comments 

As a general qualitative comment to the experimental 
patterns reported so far, it is worth to mention that the 
measured critical magnetic fields, H'^^ reported in the fig- 
ure captions, are systematically higher than those in Ta- 
ble I expected from Eq.Q (2-3 times for the Nb-Nb 



samples and 3-5 times for the Nb-NbN ones). One might 
conclude that the expression for H'^ is wrong or that 
the parameters used for the calculation are unreason- 
able. Nevertheless, when comparing the measured and 
expected values, ^q/tiC t^h^'j = ii^H'^d'^/ C t^h^ of the 
critical control currents, we found a quantitative agree- 
ment within the experimental uncertainties. In fac t, set - 
ting Kt = l\n Eq.([2Qa|) and = a = 3.3 in Eq.([20b]), 



we expect the top and bottom critical control currents 
to be 0.55 and 4.7 mA, respectively; these values have to 
be compared with those found in the experiments and 
stated in the caption of Fig. |6] We believe that the ex- 
pression for the critical magnetic field is correct, but, due 
to demagnetization effects, the externally applied mag- 
netic field, i^e, is partially screened by the electrodes^^. 



E. Flux diffraction patterns 

In a superconducting loop immersed in a magnetic 
field a current, /, circulates to expel the field from the 
loop hole. Denoting with the magnetic flux through 
the hole and with L the loop geometric (temperature- 
independent) inductance, then I = ^^/L. When the loop 
thickness, or width, are comparable with the pen- 
etration depth the magnetic flux has to be replaced by 
the London fluxoid^^ and, at the same time, the kinetic 
inductanc^^, ji^X^ /wd^ must be added to the geomet- 
ric inductance. The circulating current can be detected 
if, as shown in Fig. [7| a portion of the loop also acts 
as the electrode of an in-line LJTJ for which / is seen 
as a control current capable to modulate its critical cur- 
rent. This mechanism has been recently proposed and 
demonstrated with the aim to realize magnetic sensors 
based on LJTjJ^I^ Apart from their potential applica- 
tions, the interest for LJTJs built on a superconducting 
loop stems from the fact that they were also success- 
fully used to detect trapped flux quanta in a cosmological 
experiment aimed to study the spontaneous defect pro- 
duction during the fast quenching of a superconducting 
loop through its normal to superconducting transition 
temperatur^^. Fig. [s] shows the dependence of the junc- 
tion critical current, /c, on the external flux through 
a Nb ring of inner radius 200 /im, width 4 /im and thick- 
ness 470 nm. The magnetic flux was applied by means of 
a calibrated multi-turn coil placed underneath the chip 
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FIG. 8. Experimental flux diffraction patterns, /c($e)- The 
shielding currents induced by the external flux circulate in the 
top electrode: Ic,max = 2.25mA and $c = dzTOO^o- 



holder whose axis was perpendicular to the loop plane. 
The pattern asymmetry is very small {a < 1%) indicating 
that the boundary conditions at the junction extremeties 
are pretty the same; this is made possible by the split- 
ting of the bias currents in the two arms of the loop. 
The static properties of LJTJs with doubly connected 
electrodes have been recently investigated^^. Since the 
junction critical field only depends on the barrier param- 
eters and not on the electrode configuration, in keeping 
with our previous findings, it is easy to derive that the 
critical magnetic flux, is proportional to the loop in- 
ductance, L, and 



^0 



(30) 



same inductance per unit length. Following Weinhact we 
have generalized the conditions at the junction bound- 
aries for those more realistic cases in which the electrode 
widths are the same, but their thicknesses and materials 
are unlike. This case requires the separated knowledge 
of £5 and Ct^ the inductances per unit length of, respec- 
tively, the base and top electrode, related to the magnetic 
energy stored within a London penetration distance of 
the film inner surfaces. One interesting feature is that 
the inductance ratio is directly related to the so far un- 
explained asymmetry in the magnetic diffraction pattern 
of symmetrically biased in-line junctions. Later on the 
modeling was extended to the more common situation 
in which also the electrode widths can be different and 
useful expressions for £5 and Ct have been proposed, as 
far as the film widths are not much wider than the tun- 
nel barrier. Our approach also include junctions with 
a mixed in-line and overlap bias configuration. We like 
to stress that our analysis on windows- type junctions is 
restricted to the cases when the film widths are larger 
than, but comparable with, the junction width. We have 
reported an extensive experimental study of the static 
properties of long in-line junctions having different ma- 
terials and various geometrical configurations. One more 
interesting feature was that the junction critical current 
is equally modulated by an external magnetic field or a 
control current injected into any of the electrodes. The 
same behavior is reproduced in presence of a magnetic 
flux linked to a doubly connected electrode. We have 
shown that our experimental data are compatible with 
Eqs.(20a) and (b) that have to be taken into considera- 



tion in the assessment of the self-field effects in Josephson 
devices. 



clearly, it is assumed that the presence of the tunnel junc- 
tion does not change significantly the inductance of the 
loop. Eq.(30) can be exploited to determine the loop 



inductance in all those cases where the loop has an ir- 
regular geometry or is affected by the presence of other 
superconducting elements. We have successfully tested 
Eq.(30) for several circular and rectangular loops with 



different dimensions. In all cases a small portion of the 
loops constituted the top electrode of a LJTJ whose base 
electrode was made by a larger Nb film. As expected, no 
significant difference was found between the inductance 
values obtained for Nb and NbN loops. This made us 
confident with the reliability of the whole experimental 
procedure and with the values chosen for the bulk mag- 
netic penetrations of these two materials. 
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V. CONCLUSIONS 



We have addressed the issue of the self-field effects 
in long Josephson tunnel junctions which were tradi- 
tionally used to investigate the physics of non-linear 
phenomena^^. In all previous works on LJTJs it was 
implicitly assumed that the junction electrodes had the 
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